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ABSTRACT 


This report discusses advances in the development of a multi-frequency 
modulation (MFM) packet communications system on an industry standard 
computer. Transmitter and receiver programs are described that control vector 
signal processors and data acquisition devices. Further, these programs encode, 
modulate, demodulate, and decode the MFM signal. The throughput data rate 
was doubled, the encoding/decoding process was implemented in near real-time, 
and a personal computer plug-in board was designed and built to provide syn- 
chronization between the transmitter and receiver. This MFM implementation 
provided acceptable bit error rates at input signal-to-noise ratios of 15 dB and 
above. 
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Some of the terms used in this thesis are registered trademarks of commercial 
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of these programs without additional verification is at the risk of the user. 
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I. INTRODUCTION 


Multi-frequency modulation (MFM) provides many advantages for typical 
digital communication links. The entire system is implemented on industry 
standard personal computers (PC). Analog devices, such as voltage multipliers 
and phase-locked-loops, are not required in the modulation of the baseband sig- 
nal or in the demodulation of the received bandpass signal. The frequency range 
of the signal is limited only by the direct memory access (DMA) rate of the PC 
and, within this range, is easily changed under software control. Also, 
guardbands are essentially eliminated, allowing much more efficient use of the 
available frequency spectrum. Finally, the MFM technique as implemented here 
exclusively uses differentia] encoding schemes, eliminating the need for coherent 
detection. 

The scope of this thesis covers the improvement of an existing MFM com- 
munication system developed by LT Terry K. Gantenbein, USN [Ref. 1]. New 


concepts implemented during this project include the following: 


e the use of vector signal processors to make the signal encoding and decoding 
process near real-time, 


@ the design of software encoding algorithms to double the data rate of the 
existing System, and 


e the design and implementation of a PC plug-in board used to synchronize 
the receiver to the transmitter. 


The theory of MFM is briefly reviewed in Chapter II, the system design is 


discussed in Chapter III], the implementation is delineated in Chapter IV, and the 


analysis of the system is presented in Chapter V. The reader interested in a more 


thorough discussion of the properties of MFM is referred to Refs. 1 and 2. 


Ii. BACKGROUND OF MULTI-FREQUENCY MODULATION 


A. THE MFM SIGNAL PACKET 
The MFM signal packet has time and frequency characteristics which give it 
unique properties. The following terms are used in the description of the MFM 


signal packet: 


T : Packet length in seconds 
AT : Baud length in seconds 
k, : Baud length in number of samples 
by : Phase of the k™ tone in the /” baud 
Ay, : Amplitude of the k™ tone in the /" baud 
L: Number of bauds per packet 
At: Time between samples in seconds 
J, = 1/At: Sampling frequency or A/D and D/A conversion rate in Hz 
Af = 1/AT: Frequency spacing between tones 


K: Number of MFM tones in a baud 


These MFM parameters are depicted in Figure 1. 

The MFM signal packet is composed of a series of “bauds” of length AT . 
By choosing AT = k,At, each baud consists of k, time samples spaced Ar seconds 
apart. The vertically stacked bins in Figure 1 represent the frequencies (tones) 


in each baud. Each baud has the same tones, but the tones are encoded with 


‘different information from baud to baud. One property of MFM is that the 
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Figure |. Time/frequency parameters of MFM signal packet 


frequencies of each of the tones in a baud are integer multiples of Af, thus the 
tones are harmonics with frequency spacing Af Hz [Ref. 2: pp.6-7]. In 
Figure 1, tones f, and y, are (k, —k,)Af Hz apart. The sampling theorem requires 
that the highest frequency in the signal must not exceed one-half the sampling 


rate, f./2. Thus, in an MFM Signal packet the highest representable frequency is 


In Figure 1, symbol ‘/k’ is a single tone in the baud with frequency kAf, 
amplitude A,, and phase @,. The amplitude and phase of the tone are assigned 
in accordance with any Suitable encoding scheme. Selection of a scheme to rep- 
resent four bits per tone is a focal point of this thesis and is presented in Chapter 

I]. The frequency range of the signal is assigned by choosing the-harmonic val- 


ues k, andk, to coincide with the desired passband. Herein lies one of the 


flexibilities of MFM. Given a sampling frequency, the desired passband is 
selected by simply assigning the correct values for k, and k;. 
The mathematical description of the analog MFM signal in any baud is given 


by [Ref. 2: pp.6-7] 


k,[2-1 


x)= )) Ay cos2nkAft+ by), (I- NAT Sts IAT (1) 
k=1 


and the sampled discrete time representation is 


k,|2~1 


x;(n) = >, Ax cos( aie + Or) O<n< k, — 1. (2) 
x 
k=] 


B. MFM SIGNAL PACKET CONSTRUCTION AND MODULATION 

The Inverse Discrete Fourier Transform (IDFT) with its symmetry properties 
is ideally suited for the generation of the MFM signal packet. A symmetry 
property of the Discrete Fourier Transform (DFT) is the following: the DFT of 
a k,-point real-valued sequence is a k,-point complex sequence with conjugate 
image symmetry. That is, the values in the second half of the complex sequence 
will be the complex conjugates of their image values in the first half of the com- 
plex sequence. This relationship is diagrammed in Figure 2 for k, = 16 points. 
Conversely, the IDFT of a conjugate symmetric array of frequency domain val- 
ues produces an entirely real sequence of time domain values. This property is 


exploited as follows. After all the tones have been encoded, the amplitudes and 
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Figure 2. IDFT data structure 


phases are converted from polar to rectangular coordinate form. Next, these 
rectangular frequency plane coordinates are loaded into the first half of a k,-point 
compiex array and the complex conjugates of the values in the first half of the 
array are placed in the second half of the array at the image frequencies. The 
IDFT is performed on this array, yielding k, real time-domain values. This proc- 
ess is repeated for each baud, giving k,Z real samples. At this point, the packet 
is ready to be modulated. 

Modulating the packet is simpler than creating it. Since all the frequency in- 
formation was assigned in the selection of the passband, all that must be done to 
modulate the packet is to pass it through a digital-to-analog (D/A) converter at 
clock rate f. The block diagram of the MFM transmitter is given in Figure 3. 


The resultant D/A output contains the desired frequencies in the range initially 
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Figure 3. Block diagram of the MFM transmitter 


selected, k,Aftok,Af. The passband used in this thesis models an acoustic chan- 
nel in the range 16-20 KHz; the approximate MFM parameters are shown in 


Table 1 below. 


Table 1. PARAMETERS FOR A 1/15TH SECOND MFM SIGNAL PACKET 
IN A 16-20KHZ PASSBAND 


oven | or | vw [ow [ ve [me [os 
Ea 
Teme | ffm |@|* |= 
fesse | | @ fo [om | [om 
msc | | [| sw [aw [ 
> [= pe fo [a 


Note that any reference to the “size” of a baud describes both the number of 


samples in the baud and the number of points in the DFT used to generate the 
samples in the baud. The reader interested in voice-band channel applications is 


referred to Salsman [Ref. 3]. 


C. THE MFM RECEIVER 

The MFM reception process is the reverse of the MFM transmission process. 
After the receiver is synchronized to the transmitter (see Chapter II]), the data 
is sampled through an analog-to-digital (A/D) converter at the sampling fre- 
quency f,. A k,-point DFT is performed on k, samples. The second half of the 
array of frequency domain values is discarded as it is redundant information. 
Next, the tones are decoded and the data is presented to the information sink. 
This DF T,decoding loop is repeated ZL times until all the Lk, samples have been 


processed. The MFM receiver block diagram is shown in Figure 4. 
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Figure 4. MFM receiver block diagram 


Wl. SYSTEM DESIGN 


This chapter presents the modifications made to the MFM communication 
system developed by Gantenbein [Ref. 1]. First, the encoding process will be 
discussed with emphasis given to signal constellations, differential encoding, and 
bit rates. Next, receiver synchronization and its relationship to system phase re- 
sponse is presented. This section will also cover the development of a synchronizer 
PC plug-in board for the receiver terminal. Finally, the use of vector signal 


processors in the computation of the Fourier transforms will be reviewed. 


A. ENCODING AND DECODING 

The assignment of amplitudes and phases to the MFM tones was mentioned 
in Chapter Il. This process is known as encoding and may be based on any 
realizable signal constellation. Further, in the MFM research conducted at NPS, 
the data encoded has been files of symbols from the American Standard Code for 
Information Interchange (ASCII) character set. Though ASCII characters were 
used here, any binary information source may be used. 

Digital communications techniques generally fall into two categories, coherent 
and noncoherent. Coherent systems must extract the carrier phase information 
from the incoming signal in order to decode a received symbol, which is difficult 
and expensive. Noncoherent differential systems do not attempt to recover any 
carrier phase or amplitude information. Instead, decoding is performed by com- 
paring phase and/or amplitude differences between adjacent received symbols. 


Theoretically, MFM can be implemented using either coherent or differential 


noncoherent demodulation. For the coherent case, since the tones are all 
harmonically related, a pilot tone could be placed in a frequency bin in each baud 
several bins away from the rest of the message. The phase of this tone could then 
be recovered using traditional carrier phase recovery loops. However, this still 
may be inadequate to determine the phase of all the tones in an MFM signal in 
channels where there is significant uncompensated group delay. In the nonco- 
herent case, a differential encoding scheme is used to encode the tones; no carrier 
phase information is required. 

Though the differential schemes do not require carrier phase recovery, de- 
tection performance will be two to three dB inferior to a coherent scheme using 
the same signal constellation. For example, to achieve a probability of bit error 
of 10-7, the input signal-to-noise ratio (SNR,,) for differentially encoded phase- 
shift-keying (DPSK) is 8.5 db, while for coherent phase-shift-keying it is 6.5 db 
{Ref. 4: p. 160]. This thesis research studied noncoherent encoding/decoding 
schemes exclusively. 

The data rate (without error correction coding) in an MFM system is deter- 
mined by the number of bits encoded in each tone. Gantenbein [Ref. 1] used a 


differential quadrature-phase-shift-keyed (DQPSK) modulation constellation 


which encoded two bits into each tone. The bandwidth efficiency, B, , in any 
MFM packet communications system is given by 
No. bit 
B, = ( 9. bits ) (tones y bauds )( ] ). (3) 


baud bandwidth 


second 


For example, using the values for any of the baud sizes in Table 1, the band- 
width efficiency is two bits/sec/Hz for DQPSK which encodes two bits into each 
tone. It is apparent then that even though the MFM parameters vary, the 
bandwidth efficiency in bits/sec/Hz is the number of bits encoded in each tone. 

The challenge is to develop a Signal constellation and differential 
encoding/decoding scheme that will increase the bit rate while giving acceptable 
SNR and corresponding probability of bit error performance. To this end, two 
separate differential 16 quadrature-amplitude-modulation (DI6QAM) 
constellations are described below. 

1. D16QAM« constellation one 

The first DI6QAM signal constellation studied is shown in Figure 5. In 
this constellation, the complex plane is divided into 45° sectors and two concen- 
tric “magnitude” circles. The inner tone magnitude circle is labeled as 1X and the 
outer tone magnitude circle (which is twice the magnitude of 1X) is labeled as 
2X. The resultant complex plane has 16 distinct regions, each of which is assigned 
a four bit value. The sectors are numbered in a counterclockwise direction from 
the positive real axis starting with the inner magnitude circle. 

Using this signal constellation, the differential encoding for one symbol is 
shown in Figure 6. Assume that the previously encoded tone in a baud is at 
position A with phase 67.5° as seen in step(1). An ASCII character is retrieved 
from the message to be transmitted and the first four bits are retained. The last 
four bits are saved for the next encoding. The least significant (rightmost) bit of 
the four bit symbol is stripped off and saved, and the value of the remaining three 


bits is multiplied by 45°. The resulting value (in the example, 90°) represents the 
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Figure 5. 16QAM signal constellation 


amount of phase rotation to add to the previously encoded tone’s phase to arrive 
at the new tone’s phase. As shown here in step(3), the new tone is assigned a 
phase value of 157.5°. 

Next, the magnitude of the newly encoded tone is determined using the 
bit that was originally removed from the four bit symbol in step(2). The magni- 
tude assignment is much simpler than the phase assignment. If the least signif- 
icant bit is a zero, the magnitude of the new tone is made the same as the 
previous tone’s magnitude. If, however, the least significant bit is one, the new 
tone is assigned a magnitude different than the previous tone’s. In step (2), the 
least significant bit value is zero and the previous tone magnitude was 2X, so the 


new tone magnitude will also be 2. 


step(3) 
phase 'B’ = 157.5 degrees 


magnitude 'B’ = “large” 


step(2) 
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Figure 6. Encoding example with DI6QAM 


This constellation is referred to as the “magnitude encoding scheme” since 
the change of phase assigned to the tone is based exclusively on the magnitude of 


the three bits from each symbol. 
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In the receiver, decoding is performed to extract the tone phase and 
magnitude information. First, two adjacent tones are multiplied together using 


complex arithmetic. The phase difference A@ is found from [Ref. 1: p.11] 
ede e/Fi(¢/Fi-1)" = etig Jin — gilbi- $i). (4) 


In the absence of noise, A@ will be an integer multiple of 45°. This difference is 
then realigned to the original signal constellation by rotating it counterclockwise 
22.5°, and the bit pattern in that sector is chosen as the first three bits in the 
symbol. 

The last bit in the symbol is recovered by comparing the absolute magni- 
tudes of the two adjacent tones. Since the transmitted tones may have only one 
of two amplitudes, the decision region is placed at 1.5X, or at one-half the dis- 
tance between 1X and 2X. So if a tone’s magnitude is within one and one-half the 
magnitude of the previous tone it is decoded as the same magnitude and the last 
bit is decided to be a zero. If the tone’s magnitude is not within this range, it is 
decoded as having a different magnitude and the bit is decided to be a one. 

2. Gray-encoded DI16QAM< constellations 

Gray-encoded signa} constellations reduce the probability of bit error in 
the presence of Additive White Gaussian Noise (AWGN). In a non-Gray- 
encoded constellation, the bit patterns in adjacent sectors may differ in as many 
as three positions, as shown in sectors one and eight in Figure 5. 

The D16QAM constellation shown in Figure 7 differs from the magni- 
tude constellation in that each sector differs in only one bit position with respect 


to any adjacent sector. To encode a tone, the phase rotation to be added to the 
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Figure 7. Gray-encoded 16QAM signal constellation 


previous tone’s phase is obtained from a lookup table; there is no simple arith- 
metic relationship between the magnitude represented by the three bit pattern 


and the amount of rotation. The magnitude bit is encoded,'decoded as before. 


B. SYNCHRONIZATION 

Digital communication links require different levels of synchronization, based 
on the particular application and modulation type. MFM communication sys- 
tems require packet synchronization at the receiver. The receiver PC must have 
some method for initiating sampling at its end of the channel. This is accom- 
plished in the following manner. 

The transmitted packet contains a synchronization baud (syncbaud) ap- 


pended to its front end as shown in Figure 8. The syncbaud is a baud of size 256, 
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Figure 8. Synchronizer functional block diagram 


generated in the same manner as the other bauds in the packet. The receiver has 
a synchronization device that has as its primary feature a 128-bit correlator. The 
last 128 bits of the svncbaud sent by the transmitter are also stored in the receiver 
and used as the reference sequence for the correlator. The correlator is clocked 
at the sampling frequency J, and compares the polarities of the incoming signal 
to the reference polarities used by the transmitter. In this manner, the 
synchronizer yields a correlation peak at the last sample of the syncbaud. This 
peak triggers a latch which enables the sampling clock at the A/D converter in 
the receiver PC. Thus, sampling commences at the first sample of the first infor- 
mation bearing baud of the packet. Gantenbein [Ref. 1: pp. 24,38] implemented 
synchronization with a polarity coincidence correlator on a breadboard circuit 
external to the receiver PC. His prototype provided adequate performance but 
Was too large to implement on a PC plug-in card. 

Recent large scale integration (LSI) chip technology enabled this author to 
perform the same 128 point correlation on a plug-in card, with the added feature 


of being able 10 load the correlator reference sequence from the receiver PC. An 


overview of this device is presented here; a detailed description and schematic are 
given in Appendix A. 

The PC card-based synchronizer was developed using the TRW LSI Pro- 
ducts, Inc., TMC2221 correlator chip as its centerpiece. The TMC2221 is a 
128-point programmable digital correlator implemented in a 28 pin dual inline 
package (DIP). 

Two separate functions occur on the synchronizer PC board whose block di- 
agram is shown in Figure 9. First, the reference values are loaded into the cor- 
relation reference register on board the TMC2221 chip. Direct memory access 
(DMA) in the single byte transfer mode is used to send the reference values from 
the PC memory to the chip. The single byte transfer mode implies that a device 
requesting DMA transfer must provide a separate DMA request for each byte it 
wants placed on the PC data bus. To obtain the 128 reference values, the 
synchronizer board must make 128 byte transfer requests. The data transfer ac- 
knowledge signals returned on the PC bus by the DMA controller are then used 
to control the reference register loading in the TMC222]1 chip. 

DMaA is typically used to provide high speed data transfers from or to PC 
memory, usually on the order of thousands of bytes per second (as is done in the 
MFM transmitter D/A process to clock out up to 61440 byte samples per second 
[Ref. 1,5]). Its use here was selected with the vision that it would minimize the 
overall circuitry required to implement the transmitter and receiver hardware on 
a single PC plug-in card in the future. Further, because the reference values are 


easily changed, different acquisition delays could be tested, thus enabling 
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Figure 9. Block diagram of PC synchronizer board 


optimization of the system phase response. This concept is covered in detail in 
Chapter V, System Testing and Performance. 

The second function accomplished on the synchronizer board is the detection 
of the end of the syncbaud. After reference register loading, the correlator chip 
is ready to receive the syncbaud. The TMC2221 provides a correlation sum from 
zero to 128 on an eight-bit wide bus. By ANDing together the sum bits, any 
correlation score can be detected. This is especially important because rarely does 
the system achieve a perfect 128 point correlation. Typically, correlation sums fall 
in the range 112 to 120. This is due to system noise and channel filtering per- 


turbing the waveform so that the polarities of the received syncbaud do not cor- 
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respond exactly to the polarities of the transmitted syncbaud. Another method 
of obtaining the sampling trigger is to convert the correlator sum on the eight-bit 
sum bus through a D/A converter. The analog correlation values could then be 
compared to a threshold to determine when to start sampling the incoming 


waveform. Only the former method was implemented in this thesis project. 


C. REAL TIME SIGNAL PROCESSING 

Another goal of this thesis was to process the incoming bauds in real-time. 
For this project, real-time processing is defined as being able to complete the Fast 
Fourier Transform (FFT) and the decoding of each baud in one-half the time it 
takes to transmit it. The limiting factor here is the FFT, since the phase decoding 
is a simple table lookup function and occurs almost instantaneously in compar- 
ison with the FFT. Thus, only the FFT speed will be discussed here. 

It takes 67 milliseconds to transmit a baud size of 4096 samples when sam- 
pling at f = 61440 Hz. Using the Borland Turbo Pascal Toolbox on an IBM 
PC AT compatible machine running at eight MHz, with an INTEL 80287 nu- 
meric coprocessor, the IFFT/FFT operations each take one minute. Since up to 
eight size 4096 bauds may be received at any time in a packet, it could take up 
to eight minutes to recover the message transmitted. Clearly, this is too much time 
to pay for 32 kbytes of data. The key to making MFM competitive with high 
speed modems is to make this operation occur in real-time with vector signal 
processors (VSPs). 

Several companies have commercially available VSPs. These chips are avail- 
able as is or as the featured component in sophisticated combination data 


acquisition/signal processing boards. After investigation of several of the 
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available products, the PL1250 plug-in VSP board was purchased from 
Eighteen-Eight Laboratories. This VSP runs at 12.5 MFLOPS, using the AT&T 
DSP32 chip as its floating point engine. The board has 60 kbytes of onboard 
static RAM for storing data sent to it. The data is transferred to the board via 
either register-to-register transfers or DMA. It computes a size 4096 real FFT in 
29 milliseconds, well within the definition of real-time used here. The board is 
controlled using calls from programs written in any high level language. This card 
greatly reduced the FFT/IFFT processing time but introduced another factor 
that adversely affected the real-time processing goal. 

The data acquired and brought into the receiver PC memory had to be sent 
to the VSP. The PL1250 calls initiated this process via either register-to-register 
or DMA transfers. In register-to-register mode, the maximum data transfer oc- 
curs at 500,000 bytes per second, while DMA transfers occur at 283,000 bytes 
per second on an AT class machine at 10 MHz (Ref. 6: p.5-2]. Thus in the 
fastest of these modes, the roundtrip time for a 4096 size baud (using 32 bit real 
values to represent each sample) is 65 msecs. This is slower than the FFT com- 
putation. Future MFM research must address implementing the entire 
acquisition;signal processing algorithm on a single PC plug-in card in order to 


compcte with other high speed modem techniques. 
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IV. IMPLEMENTATION 


The architecture of the MFM packet communications system requires an ex- 
tensive software/hardware interface. Some of the required operations were writ- 
ten in Turbo Pascal, while others were coded in assembly language. It is 
worthwhile to note that the transmitter uses an NPS-designed and built D/A 
converter board that is driven with assembly language routines while the receiver 
is impiemented with the Metrabyte, Inc., DASHI6F data acquisition board and 
Pascal language calls to accomplish A/D. 

An important difference between the current transmit and receive devices is 
that the transmitter is configured to transmit up to 61440 byte samples while the 
DASHI6F can only acquire 32767 16-bit samples. The actual Metrabyte A/D 
conversion uses twelve bits, while the remaining four bits of the 16-bit word are 
used to record the A/D channel used (up to 16 channels available). This factor 
comes into play in the software because the receiver can only acquire one-half of 
the maximum length transmittable message. In this project, the message was not 
altered to accommodate the receiver; any data beyond the 32767* sample was 
ignored. 

The programs discussed here are similar in name and form to those developed 
by Gantenbein [Ref. 1]. They have been modified for the implementation of the 


VSP and to encode and decode DI6QAM. 
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A. TRANSMITTER 

The transmit program is named QAMXMIT. It uses either the magnitude 
or Gray-encoded 16QAM constellations in its encoding loop. The functions of the 
program are split into Pascal procedures (subroutines), which are shown in Fig- 
ure 10. SYNCBAUD generates the syncbaud. It uses the same parameters for 
amplitude and phase each time it is called, so the sample value polarities remain 
constant regardless of the transmitted message. Note, however, that samples of 
the syncbaud are arbitrary. These samples can be configured to establish con- 
trolled access for the message transmission and reception. 

SYNCBAUD, in turn, uses CNVTTOTIME, which calls the PL1250 
processor to perform the IFFT on the syncbaud tone phase/amplitude array. 
Next, SELECTBAUD assigns values to k,, k,, and k,, based on the desired baud 
length, such that the signal will fall in the 16-20 kHz _ passband. 
TAILORPACKET, using k,, k,, and k, from SELECTBAUD, determines the 
number of bauds that can be placed into one segment of PC memory, 64 kbytes. 
It also determines how much of the message, stored in the ASCII file 
“MESSAGE.DAT”, may be transmitted in the packet. 

DIFFENCODE performs the encoding process as described previously, in- 
cluding the placement of the complex conjugates in the second half of the IFFT 
array. It also calls CONVTTOTIME to perform the IFFT. Next, SCALEDATA 
takes the time domain samples and places them in a vector of 61440 byte values 
(the “broadcast vector”) to be clocked out through the D/A converter. 


SCALEDATA also checks each sample value to ensure that it is within the range 
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Figure 10. Block diagram of procedures in QAMXMIT 


of the eight-bit D,A converter. If the value is out of range, an error message is 
provided and the program continues to run. The encode/IFFT/scaledata loop is 
repeated until either MESSAGE.DAT is exhausted or the broadcast vector is 
filled. When the broadcast vector is ready, DMAINIT clocks the broadcast vec- 
tor through the D/A converter. DMASTOP is called to halt the DMA process 
after the 128" reference value has been placed on the PC data bus for the 


synchronizer board. 


B. RECEIVER 

The receiver program is called QAMREC. This program acquires, processes, 
and decodes a received, magnitude or Gray-encoded, MF16QAM signal. Its 
block diagram is shown in Figure 11. Exactly as in the transmitter, QAMREC 
begins by initializing the PL1250 processor. Next, GETDMABUFFER gets a 


pointer to an array which will be used to store the received samples. 
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acquire 
data i 


to 
display 


repeat if more than one baud 


Figure 11. Functional block diagram of program QAMREC 


PACKETSETUP assigns the A, and k, values and determines the maximum 
number of bauds recoverable based on ke: the receive baud length (which, of 
course, must be the same as the transmit baud length). ACQUIREDATA then 
initializes the DASHI6F A/‘D board by sending it set u, parameters via high- 
level language calls from a software package designed to drive the board (the 
software was developed by Quinn-Curtiss Software, Inc.). The parameter values 
passed to the board are the board address, the base storage address of the re- 
ceived data, DMA mode, and sampling trigger source. The DASHIG6F stores the 
received data in unipolar format, with 16 bits per sample. The four most signif- 
icant bits of each sample word are used to record the A/D channel used. 
Routine CONVERTDATA removes the DC value introduced in the 
DASHI6F sampling process and strips off the four most significant bits. The re- 


sulting 12 bit words are returned in a real-valued array and the PL1250 is called 
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to perform the FFT. The PL1250 returns values in the first half of a complex- 
valued array. DIFFDECODE retrieves the phase and magnitude difference 
between adjacent harmonics k, and k,,, and determines the ASCII representation 
for every two symbols decoded. The ASCII characters are displayed on the re- 
ceiver monitor. The CONVERTDATA, FFT, and DIFFDECODE functions 


are performed baud by baud until the entire received message is displayed. 


C. SYNCHRONIZER 

The synchronizer board is controlled by the software module, SY NCLOAD. 
SYNCLOAD calls SYNCINIT, a slightly modified version of DMAINIT. 
SYNCLOAD sets up the parameters that will be used in sending the correlation 
reference values to the TMC222]. These parameters are the number of bytes to 
transfer and the location of the bytes in memory. The real workhorse of 
SYNCLOAD is the assembly language routine SYNCINIT, which initializes the 
INTEL 8237 DMA controller on the PC motherboard to load the data from PC 
memory to the TMC2221 correlator chip. The PC bus lines driven by the DMA 
controller are used to control the loading of the reference register in the 


TMC2221. Listings of these programs are given in Appendices B and C. 
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Vv. SYSTEM TESTING AND PERFORMANCE 


The system phase response and its impact on differential phase encoding is 
described in this section. The performance measure signal-to-noise ratio out 
(SNR,,,) is defined and theoretical and observed values are compared. Last, the 
number of bit errors for 10,000 transmitted bits for various baud sizes and noise 


levels is presented. 


A. PHASE RESPONSE 

The phase responses of the system for two different cases of synchronization 
delay are presented in Figure 12. The synchronization delay is the number of 
clock cycles between the end of the syncbaud at the receiver and the start of the 
sampling. Phase response is defined as the difference between the received phase 
before decoding and the transmitted phase after encoding. It is measured for each 
tone in the baud, and represents the amount of phase shift introduced in each 
tone by the channel. It is seen here that for the 16-20 kHz passband (simulated 
by an LC filter), the magnitude of the phase difference increases with frequency 
for the zero delay case. 

If different amounts of phase change are introduced for each tone in a baud, 
there is a greater likelihood that the phase difference between adjacent tones, 
which carries the information, will be decoded incorrectly. If the amount of phase 
change added to each tone by the channel is identical, then the received adjacent 
tone symbol phase difference will only be due to the encoding process. The time 


delay, phase shift property of Fourier transforms, 
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Figure 12. Phase response for different delays 
x(t = T,) = X(fe 7", (5) 


shows that linear phase changes can be controlled by adjusting the synchroniza- 
tion delay. For a delay of one clock cycle in the sampling of the incoming 
waveform, the phase shift at harmonic number 290 is 1.8 radians. This compares 
well with the phase difference between the two curves in Figure 12, 1.9 radians 
at harmonic number 290. The discontinuity in the plot of the zero delay case is 


due to “wrap-araund” of the phase difference from —x radians to n radians. 


The phase response can be manipulated by selecting which 128 samples of the 


syncbaud are used as the correlation sequence values in the TMC2221. Note that 
the TMC2221 introduces a six clock-cycle pipeline delay between the time 
correlation occurs and the time that the correlation sum is available on the 
eight-bit correlation sum bus of the TMC222] [Ref. 7]. Thus to obtain a syn- 
chronization delay of zero it is necessary to use syncbaud values 123 through 250 
as the reference sequence. The maximum correlation sum is then output from the 
TMC2221 as the 256" sample of the syncbaud enters the TMC2221. This sum 
is then used to trigger the A/D converter to begin sampling on the first sample 
of the first information bearing baud. From measurements, however, the best 
phase response occurred when syncbaud values 124 through 251 were used as the 
reference sequence. This produces an effective synchronization delay of one. 
Thus, the phase response can be manipulated to optimize the channel for differ- 
ential encoding schemes. A synchronization delay of one was chosen for the test- 


ing described in the remainder of this chapter. 


B. SIGNAL-TO-NOISE RATIO AND BIT ERRORS 

The differentially encoded MFI6QAM system was tested for its performance 
in an AWGN environment. Different levels of input noise were introduced into 
the system and the SNR,,, was calculated for each. The SNR,,, is defined as the 
square of the normalized mean level of the complex multiplied adjacent tones di- 
vided by their variance. The decoding process involves the complex multipli- 
Cation of adjacent tones. This results in the creation of one of three tone product 
levels. The three possible product magnitudes result when adjacent tone 


amplitudes are: both small (inner magnitude), one small and one large (middle 
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magnitude), and both large (outer magnitude). The SNR,,, must be calculated 


separately for each product magnitude level. Further, though the received baud 


has the same relative tone levels (1X and 2X) as the transmitted baud, the actual 


levels are imbedded in noise. Since the actual magnitudes of these complex 


multiplications have been corrupted by noise, the means and variances have to 


be computed conditioned on the actual magnitudes of the corresponding adjacent 


transruitted tones. 


The bit error rate and SNR,,, calculations are performed in program SNR, 


which is listed in Appendix F. The general steps of the program are as follows: 


Export the magnitude and phase of the recorded encoded transmitted tones 
from the transmitter PC memory to a file in the receiver PC memory. 


For each baud, decode the received tones as previously described. 
While the received tones are being decoded, decode the transmitted tones. 


Compare the decoded transmitted and received symbols; if the received 
symbol differs from the transmitted symbol, record the number and type 
(magnitude or phase) of the bit errors. 


Combine the decoded symbols into ASCII characters and display them on 
the screen in different colors, depending on the type of error incurred in the 
symbol (if no error occurs, the character is white). 


While this joint decoding runs, sort the received complex multiplied tones 
into three bins, each representing the one of the three possible magnitudes 
resulting from multiplying adjacent tones. The sorting is done based on the 
transmitted values. 


After this sorting is completed, find the means and variances of each of the 
three bins and calculate the SNR,,,. 


The curves in Figure 13 show the SNR,,, for the inner, middle, and outer 


product magnitudes for a baud size of 1024 with various levels of SNR, . Note 


that the SNR,,, is lower as the magnitude product decreases. This is due to the 


AWGN affecting each of the product magnitude levels with equal measure, thus 
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Figure 13. SNR for 1024 baud 


causing the largest degradation when two adjacent transmitted tones are both in 
level LX. 

The SNR, curves for the inner product magnitude level are presented in 
Figure 14. Also shown on this curve is the theoretical SNR,,, for inner magni- 
tudes as given by [Ref. 8] 


2 
SNRouw = SNRi —)s (6) 
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Figure 14. Output SNR for inner powers, all baud sizes 


where c is defined as the ratio of the outer tone power to inner tone power and 
the SNR values are in linear form. In Figure 14 it is seen that the theoretical 
performance from (6) is greater than the observed values at high SNR,,. This is 
due to the system noise limiting the overall performance of the system. With no 
AWGN introduced, the maximum SNR,,, observed was about 20 dB; thus 
AWGN can only degrade the SNR,,, below this value. 

As expected, the longer bauds have generally better performance (for SNR, 
above 18 dB). This is because of the smaller frequency spacing between adjacent 


tones. At low SAR, the performance curves exceeded the theoretical. This could 


3} 


possibly be due to positive correlation of the noise of the adjacent tones, which 
would improve system performance above the theory which is based on statis- 
tically independent noise in adjacent frequency bins of the FFT. 

A more quantitative assessment of the system performance for Gray-encoded 
MFM is given in Table 2 below, which lists the number of phase and magnitude 
bit errors for 10,000 transmitted bits for each baud size at different SNR,, levels. 
The larger the baud size and the greater the SNRA,,, the better the system per- 
formance, as expected. Note that the types of bit errors are relatively evenly split 
between magnitude decoding errors and phase decoding errors. Further, because 
the phase errors are almost all of the one bit type, it is implied that most of the 
phase decoding errors were due to phases being decoded one sector away from 
their actual transmitted sectors. For this reason, the alternate (binary) bit mag- 
nitude encoding scheme was not tested, as it would clearly introduce more bit 


errors. 
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Table 2. NUMBER OF MAGNITUDE AND PHASE DECODING ERRORS IN 
10000 TRANSMITTED BITS USING GRAY-ENCODING 


TsNRINGE) [C—O CC*d 
[Mag biterrorms ‘| 8 | 4 | 0 | | 6 | 
Zit phase errors | 8 [0 | 2 [0 | 0 
3 bit phase errors | 2 [0 | 0 [0 | 0 
TSNRINGE) [SSCS SCSSC~CS~C—S—S 
[Mag bit errors | 14s] 109] 58 | 26 [8 
Tr bit phase errors | 2s] a3 oo | 19 | 26 | 
Parbit phase errors | 3 | 4 [0 | 0 | 0 | 
[3 bit phase ervors | 0 [1 [0 | 0 [0 | 
TSNRINGE) [SSC SCS 
33 [191 | 191 
ae ae 
Pabit phase errors | 0] 0 | 0 | 0 | 0 
snrin@s) | OCCCCC—S 
Tr bit phase errors | 662] 083 | 682 | 630 371 
ETRE A 
ce ee 
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VI. CONCLUSIONS AND FUTHER STUDY 


This report presents three advances which were made in a working MFM 
packet communication system. The signal processing operations were made 
real-time through the use of vector signal processors. The synchronization process 
was improved with the design and implementation of a PC plug-in card that al- 
lows experimentation with the synchronization delay and software loading of the 
reference waveform. A differential 16 quadrature-amplitude-modulation encod- 
ing technique was developed and used that doubled the data rate of the MFM 
packet communications system at NPS. The system was tested for bit error rate 
and signal-to-noise ratio performance in various levels of AWGN. Experimental 
results show that MFM in its current development can compete with other high 
speed modem techniques at input signal to noise ratios of 15 dB and higher. 

MFM has many areas in which performance can be further enhanced. 
Larger signal constellations can be developed that will increase the bandwidth 
efficiency of the MFM system. In concert with this, larger FFTs can be imple- 
mented to increase the channel data volume. The D/A conversion in the trans- 
mitter and the data acquisition process in the receiver can be expanded to 12- or 
16-bit operations to increase sample resolution which will undoubtedly be neces- 
sary for the larger signal constellations. In the receiver, to achieve real-time de- 
coding, the FFT operation should be done baud by baud as the data is being 
acquired. While a baud is being received, the FFT can be performed on the 


previously acquired baud. This is the ideal real-time decoding method and is 
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required for a continuous transmission mode. Finally, the receive and transmit 
operations should be condensed onto one PC card to enable each PC to be a 
transmit and a receive terminal. This would make it possible to implement 


networks with full duplex MFM modems. 
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APPENDIX A. SYNCHRONIZER DESIGN 


The synchronizer has two operating modes, reference register loading and 


syncbaud detection. The registers in the TMC2221 used in this project are as 


follows [Ref. 7}: 


A; - 128 point serial shift register into which the data from the hard limiter 
is placed. The data is clocked in using a 61440 Hz clock. The data in this 
register is compared with the reference sequence in R, on each clock cycle to 
determine the correlation score. 


B, - 128 point serial shift register which is used to preload the reference syn- 
chronization values. The values are clocked in by the 61440 Hz clock used 
above. 


R, - reference latch. 128 parallel load register into which the reference values 
are placed. The reference values are loaded in from register B,;. 


The TMC2221 pins used are: 


LDR - load reference control which copies the contents of serial input regis- 
ter B,; into the reference register R; for correlation. If LDR was high on the 
previous clock cycle, the contents of input register B; are copied into input 
register R; to be used as the reference sequence. 


Al - Data input from the hard limiter is clocked in here for comparison with 
the reference sequence. 


BY - port where the reference sequence is loaded into register B;. 


RE! and RE2 - are used to select either BY or BX as the reference value 
input pin. Both are tied high to select BY as the reference input pin. 


DM7 through DMO - the correlation sum is provided on this eight-bit bus 
after a six clock-cycle pipeline delay 


REFERENCE REGISTER LOADING 


The 128-point sequence of ones and zeros that make up the reference se- 


quence are stored in the ASCII file SYNCVALS.DAT in the receiver PC. The 


values are read out to the PC bus via DMA reads in the single byte transfer 
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mode. In each DMA cycle the IOW (input/output write) line is drawn low when 
the ASCII representation of each reference value is placed on the PC bus data 
lines D7-D0. The IOW line is used to drive a latch (LS374) to hold the value of 
the DO bit position of the ASCII value which is a one for an ASCII “1”, and is a 
zero for an ASCII “0”. This latched value is then loaded into the B; register of the 
TMC2221 on the next positive going edge of the 61440 Hz clock. 

The entire reference sequence must be placed into the reference register R, 
when the 128” value is read out of the file SYNCVALS.DAT. This is accom- 
plished using the DACK3 signal which is driven low just before IOW is driven 
low on each DMA read. The DACK3 signal drives the clock input to a J-K flip- 
flop which has been configured as a toggle flip-flop (J=K =1). The output of the 
flip-flop then changes state on each DMA read and is used to drive LDR. LDR 
is high on every odd-numbered DMA read and low on every even- numbered 
DMA read, so the 128 values in B,; are parallel loaded into R; on every even- 
numbered DMA read. After the 128" DMA byte transfer, DMASTOP is called 
to halt the DMA read request. At this point, the correct reference sequence has 


been placed in R;. The schematic is given in Figure 15. 


B. CORRELATION 

After reference register loading, the synchronizer board is ready to perform 
correlation on the hardlimited signal from the link. The correlation sum is pro- 
vided in the range of zero to 255 on DM7 through DMO. Synchronization only 
occurs at the beginning of packet reception, after which any size packet may be 
continuously sampled and processed (given an MFM system that operates in 


real-time as discussed in Chapter VI). The current synchronizer board must be 
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Figure 15. Synchronizer schematic 


reset after the reception of each packet by pushing a reset switch es shown in the 


schematic. Implementation of a design for continuous reception of packets is a 


future research topic. 
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APPENDIX B. SYNCLOAD 


program syncload; 
( 


Purpose: Loads the correlation reference values into the 


T™C2221 correlator chip using DMA reads in the 
single byte transfer mode. 


Inputs: The 128 correlation reference values must be in the 


Output: The syncvals (byte values) are placed on the PC bus 
Tdinidinicinidicicididinicidinidinidich dininiicinininivicininidididicioidiiaininicicicinivieicicicicicicicick 


file SYNCVALS. DAT. 


) 

uses crt; 
typ 

reference_array = array[1..128] of byte; 
var 

: integer; 

reference_values : reference_array; 

num_ref_vals : integer; 

vals : text; 

data : byte; 


°$L SYNCINIT 


(ettnes aed aet IRI I AIR IAAI ATI NIA IAI IIA IIIA IAI AISAIIIN ) 
procedure SYNCINIT(var reference_values: reference_array; 


num_ref_vals: integer); 


external; 
(Rabid tai ti cia R IIIT I IIIA IA IAIAIN ISIS ISI IA IIIISSSAIIIN ) 
begin 
assign(vals,'syncvals. dat'); 
reset(vals); 
num_ref_vals := 127; (*declare 128 values out*) 
for j := 1 to 128 do (*read each value from *) 
begin (*SYNCVALS. DAT *) 
read(vals,data); 
reference_values[ j] := data; 
writeln(reference_values|[ j] ); 
end; 
SYNCINIT( reference_values ,num_ref_vals); (* call DMA control*) 
close(vals); (* program *) 


end. 
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APPENDIX C. SYNCINIT 


code segment word public 
assume cs: code 


public syncinit 
;procedure syncinit( BLKADDRESS : XMITPOINTER; 
: BYTECOUNT : INTEGER); 


’ 

;this procedure initializes dma channel 3 and sets the 
;parameters to output the array best by passing the address 
3;of the array start on the stack. BYTECOUNT is the number 
;of bytes to transfer and is pushed on the stack by the 
;calling program (from Ref. 1, p. 58) 


dma equ 0 
dmepage equ 80h 
syncinit proc near 
push bp 
mov bp,sp suse bp to address stack 
les di,dword ptr[ bp+6];move add of best into es: di 
mov al,4bh 3;dma chan 3 single mode, read,autoinit 
out dma+11,al 
out dma+12,al sreset first/last ff 
mov ax,es 3calc high order 4 bits of buffer area 
mov cl,4 
. rol ax,cl 
push ax ;save ax for dma start addr 
and al,0fh 
out dmapage+2,al 3store in ch 3 dma page reg 
pop ax 
and al,0f0h 
add = ax, di ;get page offset 
out dmaté6,al ;output waveform buffer start addr 
mov al,ah 
out dmat6,al 
mov ax,[ bp+4] ;output dma byte count 
out dmat7,al 
mov al1,ah 
out dmat+7 ,al 
mov al,3 ;unmask ch 3 to start 
out dma+10,al 
pop bp 
ret 6 ;pop 6 bytes off stack for addr of best 


syncinit endp 
code ends 
end 
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APPENDIX D. QAMXMIT 


program QAMXMIT; 

(GARE OTT T TIT ITII IIIT ITI IIIA IAAI TIAA 
Transmits a syncbaud and message from file ‘MESSAGE. DAT’. 

The message is differentially encoded using 16-QAM. ‘MESSAGE. DAT' 

is a text file. It should already exist before using this program. 


OUTPUT. DAT is used to collect data for testing*) 
FATA A IITA RITA IAAI IAAI IA IATA IAAI IAAT IIIA AAAI AINA IAAI ASAT ) 


uses crt,plrte55; 


const 
FIRST_ELEMENT = -28929; 


type 
TNvector = array[0..4095] of single; 
TNvectorPtr = TNvector; 
BCSTARRAY = array[ FIRST_ELEMENT. .32767] of byte; 


var 

kx, 
k1,k2,I,w, 
NUMBAUDS , MAXNUMBAUDS , 
BAUDCOUNT , BYTECOUNT, 
SYMBOLCOUNT , MAXNUMCHAR , 
MESSAGESIZE,dmachn, 
n2p,bk0psz,bkipsz, 
port ,Aadd,proc,ans3 : integer; 
Badd : word; 
MAGNITUDE, 
CHARACTERS_PER_BAUD, 
PREV_TONE_MAGNITUDE , PREV_PHASE : single; 


XREAL, XIMAG : TNvectorPtr; 

INVERSE : boolean; 

TEMPBYTE , ERROR : byte; 

BCST : BCSTARRAY; 

BYTEFILE : file of byte; 

TESTFILE : text; 

ANSWER, 

NEXTCHAR : char; 

pibuf :array[0..768] of integer; 


(*$L dmainit*) 
(*S$L cmastop*) 
(wenenee Aevnce dec eewesescuceeecesesieecscaceccccecee coe scces *) 
procedure Cnvttotime; 
(*computes inverse FFT, returns values in XREAL *) 
type 
pass = array[0..8191] of single; 
passptr = pass; 
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var 
FVALUES : passptr; 
begin 


new(FVALUES); 
fillchar(FVALUES ,sizeof(FVALUES ),0); 
for i:= 0 to kx-1 do 
begin 
FVALUES [ 2*i] 
FVALUES [2*i+1] : 


XREAL [4]; 
XIMAG [4]; 


end; 

pixfto(FVALUES ,Aadd,2*kx); 
plwtxf; 
vfieee(Aadd, Aadd, 2*kx); 
cifft(Aadd,n2p); 
cereal( Aadd, Badd,kx); 
vtieee( Badd, Badd,kx); 
plwtrn; 
pixffm(Badd,XREAL ,kx); 
plwtxf; 
dispose(FVALUES); 

end; (*Cnvttotime*) 

* 


procedure SyncBaud; 

(*Process the synchronization baud and stores 
time domain sequence at the beginning of the 
area. *) 


the 256 point 
packet storage 


J, TEMP integer; 


SYNCDATA 
SYNCMAG 
syncvals 


: byte; 
single; 


: text; 


begin 
assign(syncvals, 'syncvals.dat'); 
rewrite(syncvals); 
kx: =256; 
n2p: =8; 
SYNCMAG: = MAGNITUDE; 
fillchar(XREAL ,sizeof(XREAL ),0); 
fillchar(XIMAG ,sizeof(XIMAG ),0); 


XREAL [68]:= -SYNCMAG; XIMAG [68]:= -SYNCMAG; 

XREAL [69]:= -SYNCMAG; XIMAG (69]:= -SYNCMAG; 

XREAL [70]:= -SYNCMAG; XIMAG [70]:= SYNCMAG; 

XREAL [71]:= -SYNCMAG; XIMAG [71]:= SYNCMAG; 

XREAL [72]:= SYNCMAG ; XIMAG [72]:= ~SYNCMAG ; 
XREAL [73]:= SYNCMAG ; XIMAG {73]:= SYNCMAG ; 
XREAL [74]:= -SYNCMAG ; XIMAG [74]:= SYNCMAG ; 
XREAL [75]:= SYNCMAG ; XIMAG [75]:= SYNCMAG ; 
XREAL [76]:= -SYNCMAG ; XIMAG [76]:= SYNCMAG ; 
XREAL [77]:= -SYNCMAG ; XIMAG [77]:= ~SYNCMAG ; 
XREAL [78]:= SYNCMAG ; XIMAG [78]:= -SYNCMAG ; 
XREAL (79]:= -SYNCMAG ; XIMAG [79]:= ~SYNCMAG ; 
XREAL [80]:= SYNCMAG ; XIMAG [80]:= SYNCMAG ; 
XREAL {81]:= SYNCMAG ; XIMAG [81]:= -SYNCMAG ; 
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XREAL [82]:= -SYNCMAG ; XIMAG [82]:= -SYNCMAG ; 
XREAL [83]:= SYNCMAG ; XIMAG [83]:= SYNCMAG ; 


(*complex conjugate image*) 
for J := 68 to 83 do 
begin 
XREAL [ 256-J]:= XREAL [J]; 
XIMAG [256-J]:=-XIMAG [J]; 
end; (*for J*) 


Cnvttotime; (*compute the 256 time domain values*) 


for J:=0 to 255 do (*force values to range 0-255*) 
begin (*for d/a conversion*) 
TEMP: =round( XREAL [J] + 126); 
if TEMP < 0 then 
TEMP: =0; 
SYNCDATA: =TEMP; 
BCST[ J+FIRST_ELEMENT] : =SYNCDATA; 
end; (*for J*) 
close(syncvals); 
end; (*SyncBaud*) 
[Feseresesscinse re sel sss Sees sera Srcsecessecerscescccciss = *) 
procedure SelectBaud; 
(*SelectBaud establishes kx, k1, and k2, and n2p*) 


var 
ANSWER : integer; 


begin 
kx: =0; 
repeat 
if kx < 0 then writeln('TRY AGAIN'); 
writeln('What is the length of the bauds (kx)?'); 
writeln('i.e. 256, 512, 1024, 2048, 4096'); 
readln( ANSWER); 
case ANSWER of 
256: begin 
k1:=68; k2:=83; kx:=256; n2p:=8; 
end; 
512: begin 
k1:=135; k2:=166; kx:=512; n2p:=9; 
end; 
1024: begin 
ki: =269; k2:=332; kx:=1024; n2p:=10; 
end; 
2048: begin 
k1:=537; k2:=664; kx:=2048; n2p:=11; 
end; 
4096: begin 
k1:=1073; k2:=1328; kx:=4096; n2p:=12; 
end; 
end; (*case kx*) 
if kx = 0 then kx := -1; 
until kx > 0; 


rr ——————————————————————————————————— 


end; (*SelectBaud*) 

(Be ween new neem ween eww ween ennnnennnne Wesetesadecseccvcsasces sc *) 
procedure TailorPacket; 

(*TailorPacket sets the maximum number of baud required to 


encode the message*) 


begin 
MESSAGESIZE: = filesize(BYTEFILE); 
writeln( ‘Message is ,MESSAGESIZE,' bytes. '); 
CHARACTERS_PER_BAUD: = (k2-k1)/2; 
MAXNUMCHAR: = trunc(61440. O/kx * CHARACTERS_PER_BAUD); 
if MESSAGESIZE > MAXNUMCHAR then 
begin 
writeln('Message is to large. The last ', 
MESSAGESIZE - MAXNUMCHAR, 
' characters will not be transmitted. '); 
MESSAGESIZE: =MAXNUMCHAR; 
end; 
MAXNUMBAUDS: =trunc(MESSAGESIZE / CHARACTERS_PER_BAUD); 
if frac(MESSAGESIZE / CHARACTERS_PER_BAUD) > 0.0 then 
MAXNUMBAUDS: =MAXNUMBAUDS + 1; 


repeat 
writeln; 
writeln('Enter number of ',kx,' bauds to process. Ne 
MAXNUMBAUDS,' is the maximum. '); 
readln(NUMBAUDS); 


until NUMBAUDS in [1..MAXNUMBAUDS] ; 
end; (*TailorPacket*) 
(Hanes een ew nee n ene cen eee ee nen nee eens en enneenneeenennne ~----%) 
procedure DiffEncode; 
(* DiffEncode differentially encodes the message on a tone-to-tone 
basis. BYTEFILE is read from one byte at a time. The byte 
is isolated into two 4-bit groups. Then the first three bits 
in each symbol of 4 bits are used to determine the phase shift 
between tones, and the last bit of the 4 bit symbol is to 
determine the magnitude offset . The encoded tones are 
converted to rectangular coordinates and are stored in the arrays 
XREAL and XIMAG. Bytes partially encoded are carried over into the 
next baud by global variable TEMPBYTE *) 


var 
SHORT_VECTOR , LONG_ VECTOR, PHASESHIFT, 
TONE_MAGNITUDE , TONE_PHASE , 
PREV_TONE_ PHASE , PREV_TONE_MAGNITUDE : single; 
DELTAPHI , DELTAMAG : byte; 
J : integer; 
begin 


fillchar(XREAL ,sizeof(XREAL ),0); 
fillchar(XIMAG ,sizeof(XIMAG ),0); 
LONG_VECTOR := MAGNITUDE; 
SHORT_VECTOR := LONG_VECTOR*0. 5; 
PREV_TONE_MAGNITUDE := SHORT_VECTOR; 
PREV_PHASE >= 22.5; 


XREAL [k1] 
XIMAG [k1) 


SHORT_VECTOR * cos(22. 5*pi/180. 0); 
SHORT_VECTOR * sin(22. 5*pi/180. 0); 


writeln(TESTFILE,baudcount,' ',k1,' ",PREV_TONE_MAGNITUDE, 
* ' |PREV_PHASE); 


if SYMBOLCOUNT = 0 then 
read( bytefile,TEMPBYTE); 


for J:= (kil +1) to k2 do 
begin 
SYMBOLCOUNT := SYMBOLCOUNT + 1; 
(* seperate magnitude/phase bits *) 
if frac(SYMBOLCOUNT/2) = 0.5 then 
begin 
DELTAPHI := (TEMPBYTE and SEO) shr 5; 
DELTAMAG := (TEMPBYTE and $10) shr 4; 
end; 
if frac(SYMBOLCOUNT/2) = 0.0 then 
begin 
DELTAPHI := (TEMPBYTE and $0E) shr 1; 
DELTAMAG := (TEMPBYTE and $01); 
if NOT EOF(bytefile) then 
read( bytefile,TEMPBYTE) 
else 
TEMPBYTE := $02; 
end; 


(* differentially encode the last bit in the four bit symbol *) 


if PREV_TONE_MAG. TUDE = SHORT_VECTOR then 
begin 
case DELTAMAG of 
0: TONE_MAGNITUDE := SHORT_VECTOR; 
1: TONE_MAGNITUDE := LONG_VECTOR; 
end; 
end(* previous tone short case *) 
else (* PREV_TONE_MAGNITUDE = LONG_VECTOR *) 


begin 
case DELTAMAG of 
0: TONE_MAGNITUDE := LONG_VECTOR; 
1: TONE_MAGNITUDE := SHORT_VECTOR; 
end;(* end previous tone long case *) 
end; 


(* Now use the first three bits in the symbol to determine the amount 
of phase rotation to the next encoded tone *) 
if (ans3=1) then 


begin 
PHASESHIFT := DELTAPHI * 45.0; 
end 
else 
case DELTAPHI of 
0: PHASESHIFT := 0; 
1: PHASESHIFT := 45; 
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2: PHASESHIFT := 135; 
3: PHASESHIFT := 90; 
4: PHASESHIFT := -45; 
5: PHASESHIFT := -90; 
6: PHASESHIFT := 180; 
7: PHASESHIFT := -135; 


end;(*case DELTAPHI of*) 


(* Now assign the actual phase of the tone being encoded which is a 
function of the previous phase, and the phaseshift *) 


TONE_PHASE := PREV_PHASE + PHASESHIFT; 
if TONE_PHASE >= 360.0 then 
TONE_PHASE := TONE_PHASE - 360. 0; 


(* Now convert the magnitude and phase of the tone to rectangular 
coordinates *) 


XREAL [J] := TONE_MAGNITUDE * cos(TONE_PHASE*pi/180); 
XIMAG [J] := TONE_MAGNITUDE * sin(TONE_PHASE*pi/180); 


(* Save the newly encoded tone's magnitude and phase for the next 
encoding iteration *) 

PREV_TONE_MAGNITUDE : = TONE_MAGNITUDE; 

PREV_PHASE := TONE_PHASE; 


writeln(TESTFILE,baudcount,' ',J,' ',TONE_MAGNITUDE,' ', 
TONE_PHASE); 


end; (* end of encoding process for one tone, encode next tone *) 


(* Put the complex conjugate of the encoded tones in the second half of 
the array before computing the IFFT for this baud *) 


for J:= k1 to k2 do 
begin 
XREAL [kx - JJ:= XREAL [J]; 
XIMAG [kx - J]: =-XIMAG [J]; 
end; 
end; (*DiffEncode*) 


procecure ScaleData; 

(*ScaleData converts each real value in XREAL down to a byte 
and stores the byte in the packet storage buffer, BCST. 
INDEX establishes the location in the buffer of each byte 
in the packet. *) 


var 
INDEX,J,TEMP : integer; 
DATA : byte; 
begin 


for J := 0 to kx-1 do 
begin 
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if (XREAL [J] > 127) then 
begin 
writeln('IFFT values of of range of d/a converter’); 
halt; 
end; 
TEMP := round(XREAL [J] + 126); 
if TEMP < 0 then 
TEMP := 0; 
DATA := TEMP; 
(*256 is added to INDEX to start message bauds 
after the sync baud*) 
INDEX := J+(BAUDCOUNT-1)*kx+FIRST_ELEMENT+256; 
BCST[ INDEX) := DATA; 
(* if baudcount = 1 then 
writeln(testfile,J:4,' ',round(XREAL [J])); *) 
end; (*for J*) 
end; (*ScaleData*) 


procedure Dmastop; external; 
(*Masks DMA, stopping data transfer. *) 


procedure Dmainit(var BCST: BCSTARRAY; BYTECOUNT: integer); external; 
(*Assembly language procedure used to initialize and unmask 

the IMA for data transfer. The source code must be 

converted to a OBJ file.*) 

* 


(Be enw cme ncenennecwn nance ccwnnarwanccwccsensccccrsarenscone *) 
(Bea enn cnn w een ene nen nen enn nee ne een ren newer meee nenne *) 
begin 

dmachn: =0; 


plinit(dmachn,plbuf, sizeof(plbuf)); 
plslib('C: PL1250 PLLIB. 15‘); 


proc: =1; 
port: =$0318; 
bkOpsz: =0; 


bkipsz: =1024; 

plspre(proc,port, bkOpsz,bklpsz); 
Aadd: =$04 

Badd: aS8c00, 


new( XREAL); 
new( XIMAG); 


(*contains hex values to be encoded and transmitted*) 
assign(BYTEFILE, 'MESSAGE. DAT’); 
reset(BYTEFILE); 


(*Output file of encoded symbols. Used for system testing*) 
assign(TESTFILE, 'XMITDAT. DAT’ ); 
rewrite(TESTFILE); 


repeat 

writeln(‘enter the type of encoding (l=magnitude, 2=gray)'); 
readin(ans3); 

writeln('Enter magnitude of tones 

(greater than 65, less than 1501)'); 

read1n(MAGNITUDE); 
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until MAGNITUDE > 0.0; 
writeln('Loading sync baud. '); 
SyncBaud; 

SelectBaud; 

TailorPacket; 


SYMBOLCOUNT: =0; 
TEMPBYTE: =$00; 


writeln('Number of bauds is ' 


»numbauds); 


for baudcount := 1 to numbauds do 

begin 
DiffEncode; 
writeln('Performing IFFT ' »BAUDCOUNT, : 

NUMBAUDS- BAUDCOUNT, ' left' ; 

Cnyttotime; 
ScaleData; 

end; (*for BAUDCOUNT*) 


BYTECOUNT := 256 + ela - 1; 
rite Inebyteconnt); 


repeat 
writeln('Press return to transmit'); readln; 
Dmainit( BCST,BYTECOUNT); 
repeat 
writeln('Transmit some more? (*yes or no*) 
readin( ANSWER); 
until ANSWER in ['n','N','y','Y']; 
Dmastop; 
until ANSWER in ['n','N']; 


dispose( XREAL); 

dispose(XIMAG); 

close(BYTEFILE); 

close(TESTFILE); 
end. 
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APPENDIX E. QAMREC 
program QAMREC; 
(triete 


Purpose: This program acquires the signal,stores it in a memory buffer, 
performs the FFTs, and decodes and displays the received 


symbols in the form of ASCII text. 


Inputs : The inputs are taken form the keyboard in the form of user 
responses. The inputs are: (1) The number of bauds to be 
processed, (2) the size of the bauds, (3) the type of decoding 


scheme desired (magnitude or gray). 


Outputs: The outputs are: 


File RECDAT. DAT - the real and imaginary parts of the received 


tones 
The decoded symbols, displayed on the receive CRT. 
Joiddoivddddiniciciciinidicinkideiicnvivctiiiaodidcicicickivicictaicickk deicicicick ) 
uses Graph, Crt, tp55d16,pirte55; 
(*§1-%) 
(*$R-"r) 
const Max_Buffer = 65500; 
type 
(*TYPE for real and imaginary data for FFT routing*) 

TNvector = array[0..4095}] of single; 

TNvectorPtr = TNvector;(*Pointer for FFT data array 
which allows dynamic 
allocation of memory*) 

var 

XREAL, XIMAG : TNvectorPtr; 

ERROR, TEMPBYTE : byte; 

J,1,xradd,xroadd,proc,port, 

k1,k2,kx,ANSWER,ERR_CODE, 

BAUDCOUNT , SYMBOLCOUNT ,n2p, 

NUMBAUDS , MAXNUMBAUDS , dmachn, 

bkOpsz, bkipsz, ans2 : integer; 

MAGNITUDE , PHASE : real; 

DATAVECTOR : integer; 

DMAPOINTER : pointer; 

OUTFILE, recdat : TEXT; 

pibuf : array[0..4095] of integer; 
(Ween nec cen wenn ewww nnn nee Sewers ace cee cbdssdeccesscccescee eencance ~o*) 


procedure PacketSetUp; 


(* Defines the baud parameters given the desired baud length *) 


begin 
repeat 
clrscr; 
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if kx < 0 then writeln(' TRY AGAIN'); 
writeln('Enter baud size '); 


readin( ANSWER); 
case ANSWER of 
256: begin 
kx: = 256; 
n2p: =6; 
end; 
512: begin 
kx: = 512; 
n2p: =9; 
end; 
1024: begin 
kx: =1024; 
n2p: =10; 
end; 
2048: begin 
kx: =2048; 
n2p: =11; 
end; 
4096: begin 
kx: =4096; 
n2p: =12; 
end; 


end; (*case*) 


if kx = 0 then kx := -1; 
until kx > 0; 


MAXNUMBAUDS := trunc( (MAX_BUFFER/2)/kx); 


repeat 
writeln; 
writeln('Enter number of ',kx,' bauds to process. ', 
MAXNUMBAUDS,' is the maximum. "); 
read1n(NUMBAUDS); 
until NUMBAUDS in [ 1..MAXNUMBAUDS) ; 


kl := round(kx * 67.0 / 256.0 + 1); 
k2 := round(kx * 83.0 / 256.0); 
end; cirgeketesttp®) 
(Hew ne nce wenn ewww nem eeneenne Sesceeeacnecducsecss weaccce wamecccnccena *) 


procedure AcquireData; 


(*AcquireData initializes Metrabyte DASH-16F data acquisition 
board, using TTOOLS procedure D16_int and Dl6_ainm. Data 
transfer is controlled by the DMA controller and initialized 
by Dl6_ainm and disabled by D16_dma_int —disable. TTOOLS 
procedures a@re external procedures included by ‘uses' 
tp4d16. *) 


var 
RATE: real; 
I,CNT_NUM, MODE, CYCLE, TRIGGER, 
BASE_ADR, INT_LEVEL, DMA_LEVEL, 
BOARD_NUM, CHANLO, 


OP_TYPE, STATUS, NEXT_CNT, ERR_CODE_S : integer; 


begin 
BOARD_NUM := 0; INT_LEVEL := 7; DMA_LEVEL := 1; 
BASE_ADR := $300; 


D16_init( BOARD_NUM, BASE_ADR, INT_LEVEL, DMA_LEVEL,ERR_CODE); 
if ERR_CODE <> 0 then 
D16_print_error(ERR_CODE); 


CHANLO := 0; 
CYCLE: =0; (*O-one sweep of the DMA l-autoinitialize*) 
TRIGGER: =0; (*0 - external 1 - internal*) 
CNT_NUM: =32767; (*# of samples*) 
RATE := 10000.0;(*used for internal trigger*) 
MODE := 2; (*DMA mode*) 
3 


writeln('Ready to acquire 
Dié_ainm( BOARD_NUM , CHANLO, MODE , CYCLE , TRIGGER ,CNT_NUM, 
RATE, DATAVECTOR ,ERR_CODE); 
if ERR_CODE <> 0 then 
D16_print_error(ERR_CODE); 


STATUS := 11; 


(*status indicates the progress of acquisition. When all 
samples have been acquired status=0*) 
repeat 
D16_dma_int_status(BOARD_NUM,OP_TYPE ,STATUS ,NEXT_CNT, 
ERR_CODE_S); 
if ERR_CODE_S <> 0 then 
D16_print_error(ERR_CODE_S); 


until STATUS = 0; 
writeln('Data received’); 


if ERR_CODE <> 0 then 
D16_print_error(ERR_CODE); 

D16_dma_int_disable( BOARD_NUM,ERR_CODE); 

if ERR_CODE <> 0 then 
D16_print_error(ERR_CODE); 


end; (*Acquire*) 

(Wa wencccwcnerec cece nce were w ene cececcsecececenecceses wee wencencccccecee *) 
procedure ConvertData; 

(*ConvertData seperates channel and acquired data. and converts the 
twelve bit acquired data into turbo pascal six byte real values. These 
sample values are returned in the array XREAL []*) 


var 
AD_DATA: array[0..4095] of integer; 
I,CHAN_DATA,ERR_CODE, 
SEGMENTPART , OFFSETPART: integer; 
NEWDATAVECTOR : integer; 
TEMPPOINTER : pointer; 
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begin 
fillchar(xreal ,sizeof(xreal ),0); 
fillchar(ximag ,sizeof(ximag ),0); 
SEGMENTPART: =seg(DATAVECTOR ); 
OFFSETPART: =ofs(DATAVECTOR ) + 2 * kx * (BAUDCOUNT - 1); 
TEMPPOINTER: =ptr( SEGMENTPART , OFF SETPART); 
NEWDATAVECTOR := TEMPPOINTER; 
d16_convert_data(2048,kx,NEWDATAVECTOR ,AD_DATA[ 0] , 
CHAN_DATA,0,ERR_CODE); 
if ERR_CODE <> 0 then 
D16_print_error(ERR_CODE); 
for I:= 0 to (kx - 1) do 
begin 
xreal [i] := AD_DATA[ 4]; 


end; 
end; (**ConvertData*) 
(*-oe eS ea Geena enaanoaeneanecan See 2 OSG G2 eee FFG aGeeeeaeese noe Beane *) 


procedure DiffDecode; 


(*DiffDecode differentially decodes complex frequency domain 
arrays XREAL and XIMAG. Two decoded symbols are recombined 
into a byte and transferred to the screen*) 


var 
I : integer; 
TEMPREAL, TEMPIMAG,OLDMAG,NEWMAG :single; 
BITS , PHASEBITS , MAGBIT : byte; 
TEMNPCHAR : char; 
begin 
for I:= kl to (k2-1) do 


begin 
SYMBOLCOUNT: = SYMBOLCOUNT + 1; 


(* save the current and next magnitudes for future decoding *) 


OLDMAG := sqrt(sqr(XREAL [I] ) + sqr(XIMAG [I] )); 
NEWMAG := sqrt(sqr(XREAL [I+1]) + sqr(XIMAG [I+1])); 


(* complex multiply adjacent tones to get phase differential *) 


TEMPREAL := XREAL [I] * XREAL [I+i] + 
XIMAG [I] * XIMAG [I+]; 

TEMPIMAG := XREAL [I] * XIMAG [I+1] - 
XREAL [I+1] * XIMAG [I] ; 


(* now rotate phase by 22.5 degrees to line up with constellation 
phase sectors *) 


XREAL [1] := 0.92 * TEMPREAL - 0.38 * TEMPIMAG; 
XIMAG [I] := 0.92 * TEMPIMAG + 0.38 * TEMPREAL; 


(* decode the phase difference into the first three bits of the symbol 
to be recovered *) 


if (ans2 = 1) then (*magnitude encoding scheme*) 
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begi 


end 


it 


if 


if 


if 


n 
PHASEBITS := $00; 


(XREAL [I] > 0) and (XIMAG [I] > 0) then 
if XREAL [I] > XIMAG [I] then 

PHASEBITS := $00 
else PHASEBITS := $02; 


(XREAL [I] < 0) and (XIMAG [I] > 0) then 

if abs(XREAL [I] ) > XIMAG [I] then 
PHASEBITS := $06 

else PHASEBITS := $04; 


(XREAL [I] < 0) and (XIMAG [I] < 0) then 

if abs(XREAL {I]) > abs(XIMAG [I]) then 
PHASEBITS := $08 

else PHASEBITS := $04; 


C(XREAL [I] > 0) and (XIMAG [I] < 0) then 

if XREAL [I] > abs(XIMAG [I]) then 
PHASEBITS := $0E 

else PHASEBITS := $0C; 


else (*if phase encoding scheme*) 


begin 


end; 


(* now differentially decode the magnitudes of the tones to get the 


PHASEBITS := $00; 


if 


if 


if 


if 


(XREAL {I] > 0) and (XIMAG [I] > 0) then 
if XREAL [I] > XIMAG [I] then 

PHASEBITS := $00 
else PHASEBITS := $02; 


(XREAL [I] < 0) and (XIMAG [I] > 0) then 
if abs(XREAL [I]) > XIMAG [I] then 


PHASEBITS := $04 

else PHASEBITS := $06; 
(XREAL {I} < 0) and (XIMAG [I] < 0) then 
> abs(XIMAG [I]) then 


an 

Af abs(XREAL [I] ) 
PHASEBITS := 

else PHASEBITS := 


$0Cc 


(XREAL [I] > 0) and (XIMAG [I] < 0) then 

if XREAL [I] > abs(XIMAG [I]) then 
PHASEBITS := $08 

else PHASEBITS := $0A; 


fourth and last bit in the symbol *) 


if 


(NEWMAG > 1.5*OLDMAG) or (NEWMAG < 2*OLDMAG/3) 


then MAGBIT := $01 
else MAGBIT := $00; 


(* now jam all the bits together *) 


(*£i11 TEMPBYTE with two symbols*) 
if frac(SYMBOLCOUNT / 2) = 0.5 then 
TEMPBYTE := ((PHASEBITS or MAGBIT) shi 4); 
if (frac(SYMBOLCOUNT / 2) = 0.0) then 
begin 
TEMPBYTE := (PHASEBITS or MAGBIT) or TEMPBYTE; 
TEMPCHAR := chr(TEMPBYTE); 
write(TEMPCHAR); (* put ascii character to screen*) 
TEMPBYTE: =0; 
end; (*if frac*) 
end; (*for I*) 


begin (*main body*) 


(* initialize the PL processor *) 
dmachn: =0; 
plinit(dmachn,plbuf,sizeof(plbuf)); 
plslib('c: p1850 pllib.15'); 


proc:= 1; 
port:= $0318; 
bkOpsz: =0; 


bkipsz: =1024; 
plsprc( proc, port, bkOpsz,bklpsz); 


Get DMABuf fer( MAX_BUFFER , DMAPOINTER ,ERR_CODE); 


DATAVECTOR := DMAPOINTER; (*This statement assigns a 
generic pointer to a variable of a specific pointer 
type, i.e. integer, so that the pointer can be 
passed to the dl6_ainm routine. *) 


assign(recdat,'recdat. dat'); 

rewrite(recdat); 

writeln('Which decoding scheme (l=magnitude, 2=gray)'); 
read(ans2); 


new( XREAL); 
new(XIMAG); 


ERROR := 0; 
kx: =0; 


PacketSetUp; 


SYMBOLCOUNT: =0; 
TEMPBYTE: =0; 


AcquireData; (*AcquireData samples input analog signal*) 
(* set up address in PL FPP memory that is the starting point for where 


the acquired data is sent to for FFT processing *) 
xradd: =$0400; 
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(*begin baud by baud conversion*) 
for BAUDCOUNT := 1 to NUMBAUDS do 
begin 
ConvertData; 
plxfto(xreal ,xradd,kx); 
plwtxf; 
vfieee(xradd,xradd,kx); 
plwtrn; 
rfft(xradd,n2p); 
plwtrn; 
vtieee(xradd,xradd,kx); 
plwtrn; 
plxffm(xradd,xreal ,kx); 
plwtxf; 


for j:= 0 to kx div 2 do 
begin 
xreal [j] := xreal [2*j]}; 
ximag [j] := xreal [2*j+1]; 
end; 
ximag [0} := 0; 
for i := kl to k2 do 
begin 
writeln(recdat,baudcount,' ',I, 
XREAL [I]:10:4,' ',XIMAG [I]: 10 
end; 
Dif fDecode; 
end; 


close(recdat); 
dispose( XREAL); 
dispose( XIMAG}; 
FreeDMABuffer(MAX_BUFFER , DMAPOINTER,ERR_CODE); 
if ERR_CODE <> 0 then 
P D1l6_print_error(ERR_CODE); 
end. 


24); 


APPENDIX F. SNR 
program snr; 


(itinidinidindicicinincidiciicinciiiciicciiicciiiiccinititicciiitiiiititicttiin 

Purpose: This program is used to calculate the statistics of MFM for any 
any size baud. It determines the SNR out and counts the number 
and type of bit errors in the decoding process. 


Inputs : This program requires the following data files : 
RECDAT. DAT - the received, decoded tone real and imaginary parts 
XMITDAT. DAT - the encoded, transmitted tone magnitudes & phases 


Outputs: The statistics and number and type of biterrs are stored in the 
file output.dat. In addition, the received ASCII characters are 
displayed on the screen. If the received ASCII character is 
correct, it is displayed in white. Otherwise, it is displayed 
in one of the following colors: 

yellow - if a magnitude decoding error occurred 
green - if a phase decoding error occurred 


red ~ if a phase AND magnitude decoding error occurred 
FEI ICI IIIT IIIII III III IIIT III IIIT A TIAA T AISA AAAAAAIAAIAIHIAHIIN ) 


uses crt, graph; 


var 
answer, answer2 : char; 
i,j,n,rbaud,xbaud,rtone,xtone,dtot, 
baudcount ,numbauds ,k1,k2,kx,count, 
symbolcount ,sector,b,c,d,btot,ctot, 
badbaud,numbits , badbaud2 ,bj,cj,dj : integer; 
xoldmag,roldmag,xtempreal,rtempreal, 
xnewmag, rnewmag,xtempimag,rtempimag, “ 
totphaserrs,totmagerrs,symerrs, 
smallmag,big,sml,del,meanbig,xmagbig, bigmag, 
mbig,msml,obig,osml,mmeanbig,mmeansml, 
omeanbig,omeansm]l ,mde*.,odel,msnravg,osnravg, 
meansml,varx,snrbig,sirsml,xmagsml,snravg, 
mvarx,ovarx ,msnrbig,msnrsml,osnrbig,osnrsml, 
mxmagbig,mxmagsml ,oxmagbig,oxmagsml, 
totsnr,mtotsnr,ototsnr : real; 
xbits,xphasebits ,xmagbit ,xtempbyte, 
rbits,rphasebits ,rmagbit ,rtempbyte, 
phasebitdiff ,magbitdiff,hue,pbdl, 


pbd2 ,pbd3 : byte; 
xtempchar,rtempchar : char; 
xmitdat ,recdat,output : text; 
xreal,ximag,rreal,rimag,xmag,xphase :array[1..256] of real; 
statmat :array[1..8,1..3] of 
real; 
recdata :array[{1..48,1.. 100) 
of real; 


snrin : string[ 4]; 


(Coe en a betta dae seecocecosssecsscess *) 
procedure sort; ; 


(FAIR IIIA IIIT IIA IA IAI IAAI AAAI AIA IATA AISA AAT ARITA AAS AAAI IIASA IN 
Purpose: This procedure is used to sort the complex multiplied tones into 
three different magnitude bins, inner, middle, and small. - 


Inputs : The inputs are the adjacent transmitted magnitudes xmag[i] and 
xmag[{ i+1]}. 


Outputs: The output is the global array recdata. 
toinininiciniddinindiidt kinitcatiai ) 


begin 


if (xmag[I]=smallmag) and (xmag[I+1}=smallmag) then 
begin 
statmat[ sector,1]:=statmat[ sector, 1] +1; 
b := round(statmat[ sector, 1] ); 
recdata[ (2*sector)-1,b] : =RREAL| Ij; 
recdata[ (2*sector) ,b] : =RIMAG[ I); 
end 
else if (((xmag{I]=smallmag) and (xmag[ 1+1]=bigmag)) or 
((xmag[ I] =bigmag) and (xmag[I+1]=smallmag))) then 
begin 
statmat[ sector, 2] : =statmat[ sector, 2] +1; 
b := round(statmat[sector,2] ); 
recdata[ (2*sector)~1+16,b] : =RREAL[ I]; 
recdata[ (2*sector)+16,b] : =RIMAG[ I] ; 
end 
else (* both xmag[1] and xmag[I+1] are large *) 
begin 
statmat[ sector, 3]:=statmat[ sector, 3] +1; 
b := round(statmat[ sector,3] ); 
recdata[ (2*sector)-1+32,b] : =RREAL| I] ; 
recdata[ (2*sector)+32,b] : =RIMAG[ I]; 


begin (*main body*) 
clrscr; 
assign(output, ‘output. dat’); 
rewrite(output); 
assign(xmitdat, 'xmitdat. dat'); 
reset(xmitdat); 
assign(recdat, 'recdat.dat'); 
reset(recdat); 
writeln('Enter the input snr’); 
readln(snrin); 
writeln(‘Enter the baud leugth ‘'); 
readln(kx); - 
writeln(output,'The baud length is ',kx,' and the SNRIN =',snrin); 
writeln('Enter the number of bauds to be processed’); 
readin(numbauds); 


writeln('Enter the magnitude of the xmit short tones’); 
readin(smallmag); 
bigmag := 2*smallmag; 
writeln('Throw out any bauds 7? '); 
readin( answer); 
badbaud :=0; 
bacbaud2: =0; 
if answer in ['y','Y'] then 
begin 
writeln('Which baud ?'); 
readin( badbaud); 
writeln( ‘Any others 7‘); 
readin( answer2); 
if answer2 in ['y','Y'] then 
begin 
writeln('Enter baud #'); 
readln( badbaud2); 
end; 
end; 


(* set up for the desired baud size *) 


(* 


(* 


case kx of 


256: begin 
k1: =68; k2: =83; 
end; 
512: begin 
k1: =135; k2: =166; 
end; 
1024: begin 


k1: =269; k2: =332; 
end; 
2048: begin 
k1: =537; k2: =664; 
end; 
4096: begin 
kl: =1073; k2: =1328; 
end; 
end; (*case Kx*) 


initialize overall statistical variables *) 
TOTPHASERRS :=0; 

TOTMAGERRS :=0;SYMBOLCOUNT : =0; numbits: =0; 

pbdl: =0; pbd2: =0; pbd3 :=0; bj: =0; cj: =0; dj: =0; 
totsnr: =0; mtotsnr: =0; ototsnr: =0; 


Now, count bit errors baud by baud *) 

(* read in transmit and receive values *) 

for j:= 1 to numbauds do 

begin 
del: =0; big: =0; sml: =0; meanbig: =0; meansml1: =0; btot: =0; 
mdel: =0; mbig: =0; msm1: =0; mmeanbig: =0; mmeansm1: =0; ctot: =0; 
odel: =0; obig: =0; osm1: =0; omeanbig: =0; omeansm1: =0; dtot: =0; 
fillchar(statmat,sizeof(statmat) ,0); 
fillchar(recdata,sizeof(recdata) ,0); 
for i:= 1 to k2-ki+1 do 
begin 
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readin(xmitdat,xbaud,xtone,xmag[ i] ,xphasef i] ); 
readin(recdat,rbaud,rtone,rreal[ i) ,rimag[ i] ); 


if (xbaud <> rbaud) or (xtone <> rtone) then 
begin 
writeln('RECDAT and XMITDAT do not match'); 
halt; < 
end; (*if xbaud*) 
(* convert the xmit vals to rectangular coordinates*) 
xreal[i] := xmag[ i] *cos(xphase[ i] *pi/180); 
ximag[ i) := xmag[ i] *sin(xphase[ i] *pi/180); 
end; (*for read data files*) 
writeln; 
write(j,' '); 


(* Now commence conditional decoding of the received symbols *) 
for I:= 1 to k2-k1 do 
begin 
SYMBOLCOUNT: = SYMBOLCOUNT + 1; 
(* save the current and next magnitudes for future decoding *) 


XOLDMAG := sqrt(sqr(XREAL[I]) + sqr(XIMAG[ IJ )); 
XNEWMAG := sqrt(sqr(XREAL[ I+1]) + sqr(XIMAG[ I+1] )); 
ROLDMAG := sqrt(sqr(RREAL[ I} ) + sqr(RIMAG[ I] )); 
RNEWMAG := sqrt(sqr(RREAL[ I+1] ) + sqr(RIMAG[ I+1] )); 


(* complex multiply adjacent tones to get phase differential *) 

XTEMPREAL := XREAL[I] * XREAL[I+1] + 
XIMAG[I] * XIMAG[ I+1)}; 

XTEMPIMAG := XREAL[I] * XIMAG[I+1] - 
XREAL[ I+1] * XIMAG[I] ; 

RTEMPREAL := RREAL[I] * RREAL[I+1} + 
RIMAG[ I] * RIMAG[ I+1]; 

RTEMPIMAG := RREAL[I] * RIMAG[I+1] - ‘ 
RREAL[ I+1] * RIMAG[ I]; 


(* now rotate phase by 22.5 degrees to line up with constellation 
phase sectors *) 


XREAL[ I] := 0.92 * XTEMPREAL - 0.38 * XTEMPIMAG; 
XIMAG[ I] := 0.92 * XTEMPIMAG + 0.38 * XTEMPREAL; 
RREAL[ I] := 0.92 * RTEMPREAL - 0.38 * RTEMPIMAG; 
RIMAG[ I) := 0.92 * RTEMPIMAG + 0.38 * RTEMPREAL; 


(* decode the transmit phase difference into the first three bits of the 
symbol to be recovered * 
XPHASEBITS := $00; 


if (XREAL[I} > 0) and (XIMAG[I) > 0) then 
if XREAL[I] > XIMAG[I] then 
begin 
XPHASEBITS := $00; 
sector :=1l; 
sort; 
end 
else 
begin 
XPHASEBITS := $02; 


sector :=2; 
sort; 
end; 


if (XREAL[I] < 0) and (XIMAG[I] > 0) then 
if abs(XREAL[I]) > XIMAG[I] then 
begin 
XPHASEBITS := $04; 
sector := 4; 


sort; 

end 

else 

begin 
XPHASEBITS := $06; 
sector := 3; 
sort; 

end; 


if (XREAL[I] < 0) and (XIMAG[I] < 0) then 
if abs(XREAL[I]) > abs(XIMAG[I]) then 


begin 
XPHASEBITS := $0C; 
sector := 5; 
sort; 

end 

else 

begin 
XPHASEBITS := SOE; 
sector: =6; 
sort; 


end; 
if (XREAL[{I] > 0} and (XIMAG{[I] < 0) then 
if XREAL[I] > abs(XIMAG[I]) then 


begin 
XPHASEBITS := $08; 
sector := 8; 
sort; 

end 

else 

begin 
XPHASEBITS := $0A; 
sector := 7; 
sort; 

end; 


(* decode the received phase difference into the first three bits of the 
symbol to be recovered *) 
RPHASEBITS := $00; 


if (RREAL[I] > 0) and (RIMAG[I] > 0) then 
if RREAL[I] > RIMAG[I] then 
RPHASEBITS := $00 
else RPHASEBITS := $02; 


if (RREAL[I] < 0) and (RIMAG[I] > 0) then 
if abs(RREAL{I]) > RIMAG[I] then 


RPHASEBITS := $04 
else RPHASEBITS := $06; 


if (RREAL[I] < 0) and (RIMAG[I] < 0) then 
if abs(RREAL[I] ) > abs(RIMAG[I]) then 
RPHASEBITS := $0C 
else RPHASEBITS := $0E; 


if (RREAL[I] > 0) and (RIMAG[I] < 0) then 
if RREAL[I] > abs(RIMAG[I]) then 
RPHASEBITS := $08 
else RPHASEBITS := $0A; 


(* determine the number of bit differences between the received decoded 
phasebits and the decoded transmitted phasebits *) 
PHASEBITDIFF := XPHASEBITS xor RPHASEBITS; 
if (j <> badbaud) and (j <> badbaud2) then 
begin 
case PHASEBITDIFF of 
$01: pbdl :=pbdit1; 
$02: pbdl :=pbdi1+t1; 
$04: pbdi :=pbd1+1; 
$08: pbdl :=pbd1+l; 
$03: pbd2 :=pbd2+1; 
$05: pbd2 :=pbd2+1; 
$06: pbd2 :=pbd2+1; 
$09: pbd2 :=pbd2+1; 
$0A: pbd2 :=pbd2+1; 
$0C: pbd2 :=pbd2+1; 
$07: pbd3 :=pbd3+1; 
$OB: pbd3 : =pbd3+1; 
$0D: pbd3 :=pbd3+1; 
$0E: pbd3 :=pbd3+1; 
end; (*case PHASEBITDIFF*) 


(* now count the total number of phase decoding errors *) 
TOTPHASERRS := TOTPHASERRS + PHASEBITDIFF and $01; 
TOTPHASERRS := TOTPHASERRS + (PHASEBITDIFF and $02) 
shr 1; 
TOT>HASERRS : 
shr 2; 
TOTPHASERRS : 
shr 3; 
end; 


TOTPHASERRS + (PHASEBITDIFF and $04) 


TOTPHASERRS + (PHASEBITDIFF and $08) 


(* now differentially decode the magnitudes of the tones to get the 
fourth and last bit in the symbol *) 

if (XNEWMAG > 1.5*XOLDMAG) or (XNEWMAG < 2*XOLDMAG/3) 
then XMAGBIT := $01 
else XMAGBIT := $00; 

if (RNEWMAG > 1.5*ROLDMAG) or (RNEWMAG < 2*ROLDMAG/3) 
then RMAGBIT := $01 
else RMAGBIT := $00; 

if (j <> badbaud) and (j <> badbaud2) then 

begin 

TOTMAGERRS := TOTMAGERRS + (XMAGBIT xor RMAGBIT); 
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numbits: =numbits+4; 
end; 


(* assign colors to the text that is in error *) 
if PHASEBITDIFF > 0 then 


textcolor(138); (*1. green - phase error*) 
e if RMAGBIT <> XMAGBIT then 
textcolor(142); (*yellow - mag error*) 
if (RMAGBIT <> XMAGBIT) and (PHASEBITDIFF <> 0) then 
: textcolor(140); (*1. red - dual error*) 
if (RMAGBIT = XMAGBIT) and (PHASEBITDIFF = 0) then 
textcolor(15); 


(* now put all the bits together and color the errors *) 
(*£i11 TEMPBYTE with two symbols*) 
if frac(SYMBOLCOUNT / 2) = 0.5 then 
begin 
hue := textattr; 
XTEMPBYTE := ((XPHASEBITS or XMAGBIT) shi 4); 
RTEMPBYTE := ((RPHASEBITS or RMAGBIT) shl 4); 
end;(* if frac *) 
if (frac(SYMBOLCOUNT / 2) = 0.0) then 


begin 

if (hue = 140) or (textattr = 140) then 
textcolor(140); 

if (hue = 142) and (textattr = 138) then 
textcolor(140); 


if (hue = 138) and (textattr = 142) then 
textcolor(140); 
if (hue = 142) and (textattr = 15) then 


textcolor(142); 
if (hue = 138) and (textattr = 15) then 
textcolor(138); 
‘ XTEMPBYTE := (XPHASEBITS or XMAGBIT) or 
XTEMPBYTE; 


XTEMPCHAR := chr(XTEMPBYTE); 

RTEMPBYTE := (RPHASEBITS or RMAGBIT) or 

RTEMPBYTE; 

if (RTEMPBYTE = $20) and (textattr <> 15) 
then RTEMPBYTE := $5f; 

RTEMPCHAR := chr(RTEMPBYTE); 

(*put ascii character to screen*) 


write(rtempchar); 
XTEMPBYTE: =0; 
RTEMPBYTE: =0; 
end; (*if frac*) 
end; (*for decode xmit and rec data*) 


(* now calculate the means and variances and snrout *) 
if (j <> badbaud) and (j <> badbaud2) then 
begin 
7 for i:= 1 to 8 do 
begin 
b: =round(statmat[i,1]}); 
btot: =btot+b; 
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c: =round(statmat[i,2] ); 
ctot: =ctottc; 
d: =round(statmat[ i,3] ); 
datot: =dtottd; 
if (i=1) or (i=4) or (i=5) or (i=8) then 
begin 
for count := 1 to b do 
begin 
ee ee -1,count] ) 
+ bi 
ami: #abs(cecdatel (2*1), count] )+ 
snl; 
end; 
for count := 1 to c do 
begin 
mwbig: =abs(recdata[ (2*i)- -1+16, 
count] )+mbig 
msml: sabe( recdatal (2*1)+16, 
count) )+msm1; 


for count := 1 to d do 
begin 
obig: =abs( recdata[ (2*i)-1+32, 
count] )+obig; 
osnl: sabs(reedatal (2*i)+32, 
count] )+osml; 
end; 
end 
else 
begin 
for count := 1 to b do 
begin 
big: =abs( recdata[ (2*i) ,count] 
+big; 
sm: =abs(recdatal (2*4)- -1,count] ) 
+sml; 
end; 
for count := 1 to c do 
begin 
mbig: =abs( recdata[ (2*i)+16, 
count] )+mbig; 
msnl: sabs(reedatal (2*1)- -1+16, 
count} )+msm1; 
end; 
for count := 1 to d do 
begin 
obig: sabs(recdatal (2*1)+32, 
count] )+obig 
osml: sabs(recdatal (2*4)- 1+32, 
count} )+osml; 
end; 
end; 
end; 
if btot > 1 then 
begin 
meanbig :* big/btot; 


meansml := sml/btot; 


end; 

if ctot > 
begin 
mmeanbig : 
mmeansml : 
end; 

if dtot > 
begin 
omeanbig : 
omeansn! : 
end; 


then 


mbig/ctot; 
msml/ctot; 


then 


obig/dtot; 
osml/dtot; 


for i:= 1 to 8 do 
begin 


b: =round(statmat[i,1] ); 
c: =round(statmat[i,2] ); 
d: =round(statmat[ i,3] ); 
if (i=1) or (i=4) or (1=5) or (i=8) then 
begin 
for count := 1 to b do 


begin 
del: =del+sqr(abs(recdata[ (2*i)-1,count] )- 
meanbig)+sqr(abs(recdata[ (2*i) , count] )-meansml); 
end; 
for count := 1 to c do 
begin 


mdel: =mdel+sqr( abs( recdata[ (2*i)-1+16,count] )- 
mmeanbig)+sqr(abs( recdata[ (2*i)+16,count] )-mmeansml1); 


end; 
for count := 1 to d do 
begin 


odel: =odel+sqr(abs( recdata[ (2*i)-1+32,count] )- 
omeanbig)+sqr(abs( recdata[ (2*i)+32, count] )-omeansml); 


end; 
end 
else 
begin 
for count := 1 to b do 
begin 


del: =del+sqr(abs(recdata[ (2*i) ,count] )- 
meanbig)+sqr(abs(recdata[ (2*i)-1,count] )-meansml); 


end; 
for count := 1 to c do 
begin — 


mde1: =mdel+sqr( abs( recdata[ (2*i)+16,count] )- 
mmeanbig)+sqr( abs( recdata[ (2*i)-1+16,count] )-mmeansm1); 


end; 
for count := 1 to d do 
begin 


odel: =odel+sqr(abs( recdata[ (2*1)+32, count] )- 
omeanbig)+sqr(abs( recdata[ (2*i)-1+32,count] )}-omeansml); 


end; 


end; 
end; 
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(* use a complex magnitude in statistical calculations only if it had 
more than one occurrance *) 

if (btot > 1) then 
begin 
varx: =del/(2*btot); 
snrbig: =sqr(meanbig)/varx; 
snrsml: =sqr(meansml)/varx; 
xmagbig: =meanbig/cos( 22. 5*pi/180. 0); 
xmagsml: =meansml/sin( 22. 5*pi/180. 0); 
snravg: =sqr((xmagbig+xmagsm1)/2)/varx; 
totsnr: =totsnr+10*1ln(snravg) /1n(10. 0); 
bj: =bj+1; 
end; 


if (ctot>1) then 

begin 

mvarx: =mdel/(2*ctot); 

msnrbig: =sqr(mmeanbig) /mvarx; 

msnrsml: =sqr(mmeansm1) /mvarx; 

mxmagbig: =mmeanbig/cos( 22. 5*pi/180. 0); 
mxmagsm]: =mmeansm1/sin(22. 5*pi/180. 0); 
msnravg: =sqr( (mxmagbig+mxmagsm])/2) /mvarx; 
mtotsnr: =mtotsnrt+10*ln(msnravg)/1n( 10.0); 
cj: =cj+l; 

end; 


if (dtot >1) then 
begin 
ovarx: =odel/(2*dtot); 
osnrbig: =sqr(omeanbig) /ovarx; 
osnrsm]: =sqr(omeansm])/ovarx; 
oxmagbig: =omeanbig/cos(22. 5*pi/180. 0); 
oxmagsm1: =omeansm]/sin(22. 5*pi/180. 0); 
osnravg: =sqr( (oxmagbigtoxmagsm]1) /2) /ovarx; 
ototsnr: =ototsnrt+10*ln(osnravg)/in( 10.0); 
dj: =dj+1; 
end; 

end; 


end; (*for j := 1 to numbauds*) 
writeln(output); 
writeln(output, 'The overall inner SNROUT is ',(totsnr/bj):8:3, 


writeln(output,'The overall middle SNROUT is ',(mtotsnr/cj): 8:3, 


cis a ca overall outer SNROUT is ',(ototsnr/dj): 8:3, 
writeln(output); 

writeln(output,'Total phase decoding bit errors =', 
wricela(ontput, "pba, symbo|s with one bit phase decoding 
wrételnCoutput,'(" ,pbd2," symbols with two bits phase decoding 
sa oa a symbols with three bits phase decoding 


end. 


writeln(output,'Totel magnitude decoding bit errors * “ 
TOTMAGERRS: 5: 0); 

writeln(output,' out of ',numbits,' bits transmitted’); 
close(recdat); 

close(xmitdat); 

close(output); 
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